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these  applications  involve  a  global  Navier-Stokes  systems  with  non  standard 
boundary  conditions.  The  purpose  of  this  work  is  to  prove,  using  the  classical 
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1  Introduction 


In  this  paper  we  study  the  Navier-Stokes  equations  with  the  new  boundary 
conditions  introduced  by  the  application  of  the  Domain  Decomposition  Time 
Marching  Algorithm  (  [6],  [9],  and  [10]).  These  boundary  conditions  are  of 
slip  type  ,  they  could  appear  either  by  the  application  of  the  Domain  Decom¬ 
position  Time  Marching  Algorithm  to  a  Navier-Stokes  problem  ([9]  and  [10]), 
or  to  a  Navier-Stokes/Boltzmann  coupling  ([7]  and  [9]).  In  the  latter  case, 
they  are  similar  to  the  analytical  slip  boundary  conditions  introduced  in  [5] , 
which  are  derived  from  kinetic  theory  in  order  to  replace  solving  Boltzmann 
equations  by  solving  Navier-Stokes  equations  in  the  transitional  regime. 

We  study  here  only  the  stationary  problem.  The  treatment  of  the  time 
dependent  problem  will  follow  using  the  same  ideas  developed  here  (see  [9]), 
and  the  classical  proofs  for  the  standard  boundary  conditions  (see  [l]-[3], 
and  [4]).  We  begin  by  describing  the  strong  formulation  of  the  problem  in 
the  first  paragraph,  then  we  set  the  preliminary  results  necessary  for  our 
study.  In  the  third  paragraph  we  show  the  equivalence  of  the  strong  and  the 
weak  formulations  under  some  regularity  hypothesis,  from  which  we  deduce 
the  admissible  boundary  conditions.  We  present,  then,  the  study  of  the 
stationary  problem.  Paragraph  6,  deals  with  the  uniqueness  results  when 
the  data  are  sufficiently  small.  Finally,  we  present,  in  the  last  paragraph, 
some  conclusions. 


2  Motivations 

2.1  The  general  coupling  strategy 

For  coupling  external  Navier-Stokes  equations,  with  local  Navier-Stokes  equa¬ 
tions  (dense  regimes)  or  local  Boltzmann  equations  (transitional  regimes),  we 
introduce  two  domains,  a  global  one  f 2,  a  local  one  fly  included  in  fi,  and 
an  interface  F,-  (Fig.  1  in  which  Te  denotes  Too).  The  global  solution  W  on 
f I  and  the  local  solution  Uioc  on  fly,  are  matched  by  the  following  boundary 
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Figure  1:  The  global  geometry 


conditions,  inspired  of  Schwarz  overlapping  techniques  : 

W  =  given  imposed  value  on  Too, 

n  ■  <j(W )  •  r  =  n  •  a(Ui0C )  •  r  on  the  body  T0,  (equality  of  friction  forces) 

<  g(VF)  •  n  +  n  •  <r(VF)  •  u  =  q(Uioc)  ■  n  on  T0,  (equality  of  total  heat  fluxes) 

v  •  n  =  0  on  ro, 

Uioc  —  0  on  To,  Uloc  =  W  on  the  interface  T^. 

Above,  n  ■  a  •  n  and  n  ■  a  •  r  respectively  denote  the  normal  and  the 
tangential  force  exerted  by  the  body  on  the  flow,  with  n  the  unit  normal 
vector  to  the  body  oriented  towards  its  interior. 

The  calculation  of  Uioc  and  W  satisfying  the  above  boundary  conditions 
is  then  obtained  by  the  Domain  Decomposition  Time  Marching  Algorithm, 
which  was  introduced  by  Le  Tallec  and  Tidriri  ([6],  [9]  and  [10])  and  which 
leads  to  the  following  algorithm  : 

Initialization 

1.  Guess  an  initial  distribution  of  the  conservative  variable  W  in  the 
global  domain  Cl  ; 

2.  Advance  in  time  this  distribution  by  using  the  global  Navier-Stokes 
solver  on  N\  time  steps,  with  Dirichlet  type  boundary  conditions  on  the  body 

r0; 


3.  Deduce  from  this  result  an  initial  distribution  of  the  local  variable  Uioc 
on  the  interface  T,-  and  in  the  local  domain  fiy  ; 

4.  Advance  in  time  this  distribution  by  using  the  local  solver  on  N%  time 
steps  with  Dirichlet  boundary  conditions  on  F,  and  T0. 

Iterations 

5.  From  Uioc,  compute  the  friction  forces  n  ■  cr(Uioc)  ■  r  and  heat  flux 
q(Uioc)  •  n  on  the  body  ro  ; 

6.  Advance  the  global  solution  in  time  (N\  steps)  by  using  the  global 
Navier-Stokes  solver  with  the  above  viscous  forces  as  boundary  conditions 
on  T0; 

7.  From  W,  compute  the  value  of  Uioc  on  the  interface  r,-  ; 

8.  Using  this  new  value  as  Dirichlet  boundary  conditions  on  T;,  advance 
the  local  solution  in  time  ( jV2  steps)  and  go  back  to  step  5  until  convergence 
is  reached. 

This  algorithm  completely  uncouples  the  local  and  the  global  problems 
which  can  therefore  be  solved  by  independent  solvers. 

A  parallel  version  is  also  quite  possible  although  it  is  generally  wiser  to 
use  parallel  solvers  within  steps  6  and  8. 

Remark  2.1  The  local  problem  can  be  either  Navier-Stokes  or  Boltzmann 
equations  (see  [6]- [10]). 


2.2  The  global  Navier-Stokes  problem 

The  global  domain  ft  is  discretized  using  node  centered  cells  defined  on  an 
unstructured  grid.  Then,  at  each  time  step  n  and  for  each  cell  i,  we  solve 


r  lVn+1  -Wn  _  r 

/  - - +  E  /  Fc{Wn+1)-ni 

JCi  At  JdC.ndCj 


jev(i) 


+  /  FD(Wn+1)  -m+  I  F(Wn+1)  •  m  =  -  [  F0  •  nt. 
JdCi-T  JdCiOiToo  JdCi  nr„ 


The  fluxes  Fc  and  Fd  are  computed  at  time  step  n  +  1  and  linearized,  with 
for  example  Fc  computed  by  an  Osher  approximate  Riemann  solver  [9],  and 
[10]- 
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Figure  2:  A  boundary  cell 

On  the  body  ro,  because  of  our  special  choice  of  boundary  conditions, 
the  flux  is  given  by 

0  \ 

[  F  ■  n  =  f  n{  ■  <r{Wn+1)  ■  rii 

JdCtnT0  °  ^  JdCinTc  ni  *  &{Ui0C)  *  5 

\  q(Uloc) '  ) 

where  the  aspect  of  a  boundary  cell  (7;  is  described  in  Fig.  2. 

In  other  words,  friction  forces  and  heat  flux  are  given  explicitly  as  pre¬ 
dicted  by  the  local  solver  and  the  mass  flux  is  imposed  to  zero.  Then,  in 
order  to  have  a  well-posed  problem,  at  least  in  the  incompressible  case  (see 
next  paragraphs),  the  normal  stress  (the  multiplier  of  the  zero  mass  flux  con¬ 
straint)  cannot  be  imposed  and  must  be  obtained  from  the  solution  Wn+l. 

Remark  2.2  Imposing  friction  forces  to  the  global  solution  instead  of  no 
slip  boundary  conditions  allows  to  have  an  accurate  solution  away  from  the 
boundary  layer  even  with  a  coarse  mesh  (see  [9]). 


3  Strong  formulation 

Let  ft  be  the  domain  occupied  by  the  fluid,  T  its  boundary  as  described  in 
figure  (1);  we  assume  that  ft  satisfy  : 

ft  is  a  simply  connected  bounded  domain  in  IRjV, 

T  =  5ft  is  of  class  C1 , 

(!) 

r6  u  Too  =  r  and  rh  n  =  0, 

r f,  and  Too  are  compacts  of  nonzero  measure. 
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In  the  steady  incompressible  case,  the  global  problem  in  the  coupling 
strategy  described  above  consists  in  finding  u  and  p  satisfying  the  following 
equations  : 


— i/A«  +  uVu  +  Vp  =  /  in  fi, 

(2) 

divu  =0  in  f2, 

(3) 

u  ■  n  =  0  on  T&, 

(4) 

u  =  h  on  Too, 

(5) 

n(Vu-n)T=gT  on  Tj, 

(6) 

where  n  is  the  unit  normal  to  F,  r  is  any  unit  normal  to  n. 

We  look  for  weak  solutions  of  problem  (2)-(6),  i.e.  (2),  (3)  are  satisfied 
in  a  distribution  sense  and  equations  (4)-(6)  are  satisfied  in  some  Sobolev 
spaces. 

Remark  3.1  From  the  remark  (2.1)  we  observe  that  gT  can  be  issued  from 
a  local  kinetic  model  (Boltzmann  equations)  or  a  local  continuous  model 
(Navier-Stokes  system). 

4  Preliminary  results 

Let  Hm~2(T)  be  a  space  of  trace  functions  in  Hm(Q)  with  the  following 
norm  : 

ll’/’L-J.r  =  M{||u|U/“  €  =  4>}.  (7) 

and  V  the  closed  subspace  of  [Hx{Vt))N  defined  by 

V  =  {v  €  (H* (£l))N \divv  =  0,  v  •  n  =  0  on  Tb,v  =  0  on  Too},  (8) 
and  define 

IMI  =  llVullo,2,n.  (9) 

We  have  the  following  classical  lemma. 
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Lemma  4.1  ||.|[  is  a  norm  ofV  equivalent  to  the  norm  ||.||i,2,n- 
Lemma  4.2  ([11])  Under  hypothesis  (1),  we  have 


n  : 


r  -^ir^ 


x  n(x ) 

is  C 1  and  hence,  the  trace  operator 

(H'ityf  H2(T) 


(10) 


(11) 


u 


u  *  n 


is  continuous. 


Lemma  4,3  ([ If ])•  Let  Cl  satisfy  ( 1 ).  We  have  : 
The  injection 

whp(ci)  c  Lqi(n) 

is  compact  for  q\  satisfying  : 

1  <  qi  <  00  if  p  >  N 


1  <  qi  <  q  with  q  given  by  —  =  - 


if  1  <  p  <  N. 


(12) 

(13) 


In  particular  we  have  for  p  =  2  : 

H\tt)cLq(n)  for  2  <q<^L, 

the  injection  being  continuous.  For  N  =  2,  this  injection  takes  place  for  any 
finite  q  such  that  q  >  2. 


Lemma  4.4  ([12])  Let  Lt  be  a  bounded  connected  domain.  Then 

(i)  The  operator  grad  is  an  isomorphism  from  Lq(SI)  into  V°, 

(ii)  the  operator  div  is  an  isomorphism  from  (V0)1  into  LfifT) 
where  Lq(Q,)  and  V°  are  defined  by 

V°  =  {h€  H-\n)N,<  h,v  >=  0,  Vu  6  V}. 

L20(Ll)  =  {he  L20{ n),  j  h(x)dx  -  0}. 
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We  now  set 


.  ^  r  dvi 

v,w)=  >  /  Uj  Widx,u,v,w€V 

i,j-l  ''ft  &xj 

(14) 

<  /,  v  >=  /  fvdx,v  €  V 

Jq 

(15) 

<  9t,v  >rb=  /  grT  ■  vd'j,  v  eV 

(16) 

JTb 

((u,  u))=  /  VuVvdx,u,  v  €  V 

Jn 

(17) 

(u,  v)=  uvdx,  u,  v  €  V. 

Jci 

(18) 

If  /  is  in  (/f_1(f2));v  and  gT  in  H~2(Tb)  all  the  above  operators  (14)-(18) 
are  well  posed.  Moreover,  the  operator  defined  by  (16)  and  the  trilinear  form 
defined  by  (14)  satisfy  the  following  properties  : 

Lemma  4.5  The  mapping  from  V  into  IR  defined  by  (16)  is  continuous. 

Proof:  The  proof  follows  directly  from  lemma  4.2. 

Lemma  4.6  (i)  The  trilinear  form  b  is  continuous  on  V2  x  (V  C\  (LN( £l))N 
(ii)  Let  u,v,w  G  V.  Then  we  have  the  following  relations  : 

b(u,v,w)  =  —b(u,w,v),  (19) 

b(u,v,v)  =  0.  (20) 

(Hi)  Let  (un}n>0  be  a  sequence  in  V  such  that  un  — *■  u  weakly  in  V .  Then  : 

J.iin  i!»(un,  un,v)  =  fe('u,u,  v),Wv  €  V.  (21) 

Proof 

We  present  only  the  proof  of  (i)  and  (ii).  For  the  proof  of  (iii),  see  [14]  or 

[13]- 

(i)  Let  u  €  V  then  U{  €  Hl(Tt)  ;  from  lemma  4.3,  we  deduce  : 


which  implies 


dvi 
dx  k 

<  C'||uj||v||Vi||v||«>i||v. 


<  ||«illo,9,o||-ft—  ||o,9,h|ki||o,N,n 

OXk 


(22) 


Therefore  we  conclude  to  our  result. 

(ii)  By  definition 

b(u,v,w)  =  /  UiVjjWj. 
After  integration  by  parts,  we  obtain 


u  ■  nv  •  w 


b(u,  V,  w)  =  -  J  Vj(uiWj)ti  +  Jd 

=  -  VjUiWji  —  /  VjWjUii  +  u  ■ 

Jo  Jn  JdU 


nv  •  w. 


Since  u  6  V,  we  have  u-n  =  0  on  dCl  and  =  divu  =  0  on  Vt.  Therefore, 
we  get 

b(u,v,w)  =  —  J  VjitiWj'i  =  —b(u,w,v). 

Remark  4.1  For  N  <  4,  it  follows  from  lemma  f.3  that  V  H  (LN(Cl))N  =  V 


5  Existence  Result 

We  now  go  back  to  our  initial  formulation: 
Find  u  and  p  such  that  : 


—vAu  +  u  •  Vu  +  Vp  =  / 

in  0, 

(23) 

divu  -  0 

in  fi, 

(24) 

II 

c 

3 

on  r&, 

(25) 

XL  —  h 

on  r^, 

(26) 

< 

£ 

II 

A 

on  Tb ■ 

(27) 
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5.1  Hypothesis  on  <f>  and  h 

The  data  <p  on  Tj,  and  h  on  Too  must  verify  some  regularity  and  compatibility 


properties.  More  precisely  we  assume  that  : 

(28) 

<t>  £  Hi( iy, 

(29) 

8 

II 

CD 

(30) 

— » 

Cr 

II 

o 

(31) 

Lemma  5.1  The  relations  ( 28)-(31 )  imply  the  existence  of  u 
such  that  : 

£  (H'{ Sl)f 

V.u  =  0  in  ft, 
u  •  n  =  <p  on  Tt, 
u  =  h  on  Too. 


Proof 

Let  xp  be  a  function  defined  on  T  such  that  : 

{4>  •  n  on  Tb, 
h  on  Too. 

From  the  assumed  hypothesis,  on  F,  o  and  h  we  get  : 

4>  £  (Hi(T))N. 


(32) 

(33) 

(34) 


(35) 


Therefore,  from  the  trace  theorem,  we  can  take  u  €  //1(ft)jV  such  that  : 
u|r  =  xp.  By  using  the  relations  (30)  and  (31)  and  the  Green  Formula,  we 
obtain 

I  V  •  u  =  0. 

Jn 
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Using  the  regularity  of  u,  we  have  V  •  u  €  Lq(Q,). 

Lemma  4.4  shows  that  there  exists  v  £  Hq(CI)n  such  that  V  •  v  =  V  -  u. 
Finally,  by  taking  now  u  —  u  —  v,  we  can  see  easily  that  u  verify  the  relations 
(32)-(34). 

Now,  we  are  going  to  pick  out  some  function  £  in  Hl(fl)  such  that  : 

u  —  rot£.  (36) 


To  achieve  this,  we  will  use  the  following  proposition  : 

Proposition  5.1  ([13])  Under  the  regularity  hypothesis  (1)  and  for  N  —  3 
we  have 

rotH1^)  =  {u  £  L2(£l)N,divu  =  0,  f  u-n=t  u-n  =  0}  (37) 


Remark  5.1  This  result  is  also  valid  for  N  =  2  ([13]). 


Under  the  hypothesis  (30)-(31)  and  the  relations  (32)-(34)  we  have  : 


(38) 

(39) 


Consequently,  from  Proposition  5.1,  there  exists  £  £  i/1(U),  such  that  we 
have  (36).  From  now  on  we  assume  that  N  <  3. 


5.2  Weak  formulation 


Let  u  €  V.  After  multiplication  of  (23)  by  v,  integration  by  parts  and  by 
taking  into  account  the  boundary  conditions  (25)- (27)  and  the  relation  (24) 
we  obtain  the  following  variational  formulation  : 


Find  u0  €  V  such  that 

v((u0,  v))  +  b(u,  u0,  v )  +  h(u0 ,  u,  v)  +  b(u0,  u0,  v ) 

=  -!/(({£,  v))  -  b(u,u,v)  +  /  gTr-v+  fv, 

j  jq 

.  Vu  €  v  n  LN(n)N 


(40) 
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5.3  Equivalence  between  the  two  formulations 

We  will  show  here  that  the  variational  formulation  (40)  is  equivalent  to  the 
strong  formulation.  This  is  the  goal  of  the  following  proposition  : 

Proposition  5.2  If  ( u,p )  is  a  smooth  solution  of  (2)-(6),  then  it  satisfies 
(40).  Conversely,  ifu0  is  solution  of  (40),  then  there  exists  a  unique  function 
(up  to  a  constant)  p  £  L2{fl)  such  that  u0  +  u  and  p  satisfy  : 

(i)  (2)  and  (3)  in  a  distribution  sense, 

(ii)  (4)  in  HH Tb), 

(iii)  (5)  in  Hh( Too), 

(iv)  (6)  in  iH(T6). 

Proof 

The  direct  theorem  was  showed  in  our  introduction  of  the  weak  formula¬ 
tion.  We  then  show  only  the  converse. 

Let 

V  =  {</>  £  (Hq(Q,))n / divij)  =  0). 

Let  u0  be  a  solution  of  problem  (40)  then  after  multiplication  by  a  test 
function  in  T>( Q),  and  integration  by  parts  and  density,  we  get 

<  —  i/Au  +  (u  •  V)u  —  /,  >=  0,  €  V  (41) 

where  <,>  denotes  now  the  duality  between  the  spaces  V(Cl)N  and 
V(0,)N .  We  consider  the  operator  —grad(=  —V)  £  C{Lfy(Vt),  /f_1(fl)JV), 
and  Tt(-grad)  its  image  in  H~l{ tt)N .  We  have  the  following  lemma  (  see 
[12]) 

Lemma  5.2  7 Z(-grad)  is  identical  to  V°  with  V°  defined  by  : 

v°  =  {he  H-'ittf/  <  h,i/>  >=  0,V</>  €  V}. 
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The  application  of  this  lemma  implies  the  existence  of  p  €  To(fi)  such 
that 

—uAu  +  u  ■  Vu  —  /  +  Vp  =  0  in  (/T^fl))*. 

As  u  6  V,  u  ■  Vu  and  /  6  I/2,  then 

—i/A u  —  Vp  =  div{— i/Vu  +  p/d) 


=  (/-u-Vu)<e(L2(12))n 
Applying  the  Green  formula  implies  that  there  exists 


q  =  (— i/Vu  +  pld)  ■  n  <E  H  2  (512) 

such  that 

/  —  div(— z/Vu  +  pld)  ■  v 
Ju 

+  /  z/VuVu  —  /  pdivv  —  [  q  ■  v,  Vv  £  (H1(  12))^. 

Jn  in  ./an 

By  replacing  diu(— i/Vu  +  pld)  by  (/  —  u  •  Vu)  we  get 

i/((u,  u))  +  6(u,  u,v)  —  J  pdivv  -  q  ■  v  -  J^fv  =  0 


Vu  €  (ff1^))* 


By  subtracting  (43)  from  (40),  it  remains 


/  9  *  v  “  L 

/an  4rfc 


/  ^t(v-t) 


=  0,  Vu  €  V. 


Now  let  zu  arbitrary  element  of  //020(r6),  where  //^(Tfc)  is  defined  by  : 


^(Tfc)  =  {u  €  T2(T6);  3zu  €  H1^),  zu  =  0  on  Too  andzu  =  u  on  Tfc} 
(see  [11]  for  more  details  about  this  space). 
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Let  v  G  V  such  that 


v-n  =  0  on  r6, 
v  ■  t  =  w  on  r6, 


q-T  =  t  •  (— i/Vu  +  pld )  •  n 
=  gT  in 

Remark  5.2  We  cannot  impose  the  normal  stress  in  (3)-(7)  because  of  the 
presence  of  a  pressure  term  in  q. 

5.4  Existence  result 

Now,  we  are  able  to  set  out  the  main  result  of  this  section. 

Theorem  5.1  Let  f  G  H~1(Cl)  ;  then  there  exists  u0  G  V  solution  of  (40). 

Proof 

We  proceed  as  in  the  usual  proof  of  existence  for  the  incompressible 
Navier-Stokes  equations.  Let  us  define 

[f>m(u),u]  =  v\\u\\2  +  z/((u,u))  +  b(u,u,u)  -  I  grTU 

—  fu  +  b(u ,  u,  u). 

J  Q 

By  using  the  continuity  of  b  we  obtain  : 
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(43) 


b(u,u,u)  <  ci||u||0,g,n||'S||i,2,fi||«||o)Ar 
with  q  satisfying  (12)-(13).  In  addition  we  have  : 

b{u,u,u)  <  ca||u||?j2>0||«||.  (44) 

Using  lemma  4.1  and  lemma  4.2  we  arrive  to  : 

[V>m(u),u]  >  v\\u\\2  +  b(u,U,u) 

HMIMNI  +  c2||«lll,2,n  +  +  ||/||-l,2,n)- 

So,  to  be  able  to  use  the  Brouwer  fixed  point  theorem,  it  is  sufficient  to 
show  that  we  have  for  some  /?  >  0,  the  following  relation  : 

u\\u\\2  +  b(u,u,u )  >  /?||u||2,Vu  £  V.  (45) 

This  is  the  goal  of  the  following  lemma  : 

Lemma  5.3  V  7  >  0  we  can  choose  u  verifying  (32)-(33)  such  that 

|6(u,«,u)|  <  7IMI2-  (46) 

From  the  above  lemma,  we  conclude  to  the  existence  of  a  given  k  such 
that 

hM«),(u)]  >  0, 

for  any  u  6  V  with  ||u||  =  k.  Then  the  Brouwer  fixed  point  theorem  can  be 
applied  to  the  function  tpm(u)  inside  any  ball  of  radius  k  belonging  to  any 
given  finite  dimensional  approximation  Vm  of  V .  From  this,  we  deduce  the 
existence  of  a  Galerkin  approximation  «m  €  Fm  of  the  solution  u  of  (40).  By 
standard  compactness  arguments,  it  follows  that  um  weakly  converges  to  a 
solution  u0  of  (40)  in  V.  (see  [5]  for  more  details). 
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6  Uniqueness  results 


Under  the  assumption  that  iV  <  3,  we  have  the  following  result  : 

Theorem  6.1  Assume  that  u  in  (ZN(Si))N  is  sufficiently  small  so  that 

\b(v,u,v)\  <  ||M|2,Vu<=  V,  (47) 

and  v  is  sufficiently  large  so  that 


^  >  m\f\\v  +  \\9r\UTb) 


where  C  is  the  constant  in  (22)  and 


/  =  vAu  -  EjUj—u  +  / 


and 


(48) 


d  . 

9r  =  +  9r 

Proof 

We  proceed  as  in  the  proof  for  the  standard  homogeneous  boundary  con¬ 
ditions  (see  for  instance  [13]). 

Suppose  that  u0  is  solution  of  the  problem  (40).  Taking  v  =  ua  in  (40) 
and  using  (20),  we  obtain  : 


v\\un 


\v  =  ~Kuoi  u,  u0 )  -  i /((«,  u0))  -  h(u ,  u,  u0 )  +  jf  grT-u0  +  fu0,  (49) 


Using  the  Green  formula,  we  obtain 


v\\u0 


2 

V 


d 

-b(u0,  u,  u0)+  <  vAu  -  EjUj—u  +  f,u0> 


+ 


d  . 

<  v-^u  +  gT,u0  >rb 


(50) 

(51) 


=  -b(u0,  u,  u0)+  <  /,  u0  >  +  <  gr,  u0  >r6  (52) 
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Since  /  €  and  u  G  (tf1^))^,  An  G  {H-\tl))N  and  then 

/  G  (tf_1( fl))*.  Similarly,  since  gT  G  (tf~5(r6))^  and  u  G  (i71(H))Ar  we 

have  —  G  (//^(Tj))^  and  then  <7T  G  (tf-i(Pj))w,  and  the  equality  above 

(JTl 

is  well  defined.  The  rest  of  the  proof  is  an  adaptation  of  the  proof  for  the 
standard  homogeneous  boundary  conditions  (see  for  instance  [13]). 


7  Conclusion 

In  this  paper  we  have  shown  that  the  slip  boundary  conditions  resulting 
from  the  application  of  the  Domain  Decomposition  Time  Marching  Algo¬ 
rithm  to  either  Navier-Stokes/Navier-Stokes  coupling  ([9]  and  [10])  or  Navier- 
Stokes /Boltzmann  coupling  ([9],  [7],  and  [8])  are  admissible  and  lead  to  a  well 
posed  global  problem  at  least  in  the  incompressible  case. 
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